Measurement of neutron yield by 62 MeV proton beam on a thick Beryllium target 



M. Osipenko", M. Ripani", R. Alba^ G. Ricco'^, M. Barbagallo^ P. Boccaccio'^, A. Celentano% N. Colonna^ 
L. Cosentino^, A. Del Zoppo^, A. Di Pietro'', J. Esposito'*, P. Figuera'', P. Finocchiaro'', A. Kostyukov-'', C. Maiolino'', 

D. Santonocito^ M. Schillaci^ V. Scuderi^ CM. Viberti^ 

" INFN, sezione di Geneva, 161^6 Geneva, Italy, 
' INFN, Laboratori Nazionali del Sud, 95123 Catania, Italy, 

" INFN, sezione di Ban, 70126 Ban, Italy, 
INFN, Laboratori Nazionali di Legnaro, 35020 Legnaro, Italy, 
Dipartimento di Fisica dell'Universita di Genova, 16146 Genova, Italy, 
^ Moscow State University, Moscow 119992, Russia 



m 
o 

00 
(N 



X 



> 

(N 

(N 
O 
m 



X 
S3 



Abstract 

The design of a low-power prototype of neutron amplifier recently proposed within the INFN-E project indicated the 
need for more accurate called for detailed data on the neutron yield produced by a proton beam with energy of about 70 
]V[eV impinging on a thick Beryllium target. Such measurement was performed at the LNS superconducting cyclotron, 
covering a wide angular range from to 150 degrees and a complete neutron energy interval from thermal to beam energy. 
iNeutrons with energy above 0.5 MeV were measured by liquid scintillators exploiting their Time of Flight to determine 
the kinetic energy. For lower energy neutrons, down to thermal energy, a '^He detector was used. The obtained data 
,are in good agreement with previous measurements at degree with 66 MeV proton beam, covering neutron energies 
'>10 MeV, as well as with measurements at few selected angles with protons of 46, 55 and 113 MeV energy. The present 
iresults extend the neutron yield data in the 60-70 MeV beam energy range. A comparison of measured yields to MCNP 
and Geant4 Monte Carlo simulations was performed. 

Keywords: neutron yield, thick target 
•PACS: 29.25.Dz, 29.30.Hs 



1. Introduction 

Research on fast neutron reactors addresses, among 
others, the problems of the limited availability of fissile ma- 
terials and the difficulties in nuclear waste disposal. This 
is because fast neutrons induce fission on larger number 
of nuclides, including 238-(j g^^j Minor Actinides like e.g. 
^'^-'^Am, while they undergo much less capture processes 
leading to the build up of high toxicity waste. Accelerator 
■Driven System (ADS) [l| represents a promising solution 
to this problem. It consists of a fast sub-critical reactor 
fed by an external neutron source of sufficiently high in- 
tensity, that can be generated by a particle accelerator or 
by a fusion reactor. 

The high power proton cyclotron that will be installed 
at Laboratori Nazionali di Legnaro of INFN for the SPES 
project [H offers the possibility to build a low-power ADS 
prototype for research purposes. The design of such proto- 
type has been developed by a broad collaboration between 
INFN and other bodies and is documented in Ref. [s*]. It 
is based on a 0.5 mA, 70 MeV proton beam impinging 
on a thick Beryllium target and a fast sub-critical core 
consisting of solid lead matrix and 60 UO2 fuel elements, 
enriched to 20% with ^-^^U. Both, production target and 



reactor core, are cooled by continuous helium gas flow. 
The ADS is expected to have the effective neutron multi- 
plication factor kf.ff = 0.946, neutron flux (/) = 3-^6xl0^^ 
n/cm^/s and thermal power of 130 kW at 200° C. Such a 
facility would allow to study the kinetics and dynamics of 
the fast reactor core, the burn-out and transmutation of 
radioactive waste, as well as issues related to system safety 
and licensing. 

The choice of the production target material is deter- 
mined by the low beam energy of the cyclotron. ^Be bom- 
barded by protons, in fact, provides an abundant neutron 
source whose spectrum is sufficiently hard to burn minor 
actinides. The design of the proposed ADS requires careful 
knowledge of the neutron yield produced by the 70 MeV 
proton beam on a thick ^Be target. The existing data in 
the given energy range are rather scarce and not compre- 
hensive. The integrated yield at 70 MeV was measured in 
Ref. Q, while differential yields for few angles were mea- 
sured at various beam energies in Refs. (5_, 6, 7, 8, 9, lQi,[ll|. 
However, only the data from Ref. have sufficiently sim- 
ilar beam energy, also the results are limited to de- 
grees and neutron energies >10 MeV. This lack of data 
demanded a dedicated measurement of neutron yield pro- 
duced by a proton beam on a thick ^Be target. The mea- 
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surement was performed at the Laboratori Nazionali del 
Sud (LNS) of INFN using the existing superconducting cy- 
clotron [13) ■ The cyclotron was usually operated with 
a 62 MeV proton beam, hence we also selected this, well 
tuned, energy, close enough to the project goal. Time of 
Flight in liquid scintillators was used to measure neutrons 
with energies >0.5 MeV, while a ^He detector has been 
used for neutrons of lower energy. The precision of the ex- 
periment was of the order of 10%, dominated by system- 
atic uncertainties on the detector efficiencies and absolute 
normalization. 

In the first three sections we discuss experimental setup, 
detector geometry and energy calibrations. In further sec- 
tions the data analysis is presented. Finally in the last few 
sections the obtained data are discussed and compared to 
previous data and Monte Carlo simulations. 



made of two Teflon wheels 2 mm thick with outer diameter 
equal to the beamline pipe size and inner diameter fitting 
the target cylinder as shown in Fig. [TJ The front face of the 
target with respect to the beam had a guard ring polarized 
at -200 V to prevent (5-electron emission from the target. 
Tests made at different values of bias voltage showed no 
significant variation of the collected charge. 




2. Experimental Setup 

The 62 MeV proton beam from the LNS superconduct- 
ing cyclotron was delivered in the experimental hall with 
a beam current of 30-50 pA. The cyclotron beam struc- 
ture (RF=40 MHz) was modified by suppression of four 
bunches out of five. This way 1.5 ns wide beam bunches 
arrived on the target with period of 125 ns. The target 
consisted of a solid 3 cm thick ^Be cylinder with 3.5 cm 
diameter. The target thickness and width were chosen 
to ensure complete absorption of the proton beam. The 
electric charge deposited by the beam on the target was 
measured by a digital current integrator and used for ab- 
solute normalization of the data. Two different current 
integrators were used during the experiment for compari- 
son. 

The precision of the present experiment depended on 
careful reduction of environmental background. To reduce 
7 production and rescattering of the large fiux of pro- 
duced neutrons from surrounding materials specific mea- 
sures were taken. First of all the beam line was located 
at a height of 168 cm above an empty fioor. All beamline 
components downstream of the target were completely re- 
moved. Closest wall of the experimental hall was located 
at more than 4 m distance from the target. During prelim- 
inary measurements, the 3 mm thick steel beamline proved 
to be a significant source of background, in particular of n- 
7 conversion. This was also confirmed by MCNP fl3| sim- 
ulations. To reduce this contamination a new, lightweight, 
final segment of the beamline was installed. It consisted of 
100 cm long carbon fiber tube with 3 mm thick walls man- 
ufactured by Quick Batten [13]. The tube was connected 
to the main beamline via an ISO-K Aluminum fiange glued 
to the carbon fiber. At the downstream side the tube was 
sealed by a glued 2 mm thick Aluminum endcap. Such so- 
lution permitted to reduce the background by more than 
a factor of 2, maintaining the vacuum quality better than 
10^^ mbar. The target holder was also specially designed 
to avoid large background production and good electric in- 
sulation to guarantee the best charge collection. This was 



Figure 1: (color online) Target holder drawing. The beryllium target 
is represented by the yellow cylinder. 



The neutrons produced in the target were measured by 
the ToF technique. To this end, the RF signal from the 
cyclotron was used as the reference time. Eight neutron 
detectors were installed simultaneously around the target, 
at the same height of the beamline at different angles and 
at two different distances (150 cm and 75 cm) around the 
target as shown in Fig. [51 Many separate runs were taken 
during the experiment. For each run all eight detectors 
were moved simultaneously to other angles for a total of 
16 positions, sometimes repeated for consistency checks. 
Near (small) detectors, installed at half distance (75 cm) 
from the target with respect to far (large) detectors, were 
used to measure low energy neutrons (r„ <2 MeV). 
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Figure 2: Schematic view of the experimental setup. The detec- 
tors were moved simultaneously at other angles to cover 16 angular 
settings from to 150 degrees. 



All detectors were installed on light Aluminum tubes 
by using Tefion holders. The tubes were mounted on an 
Aluminum foot for mechanical stability. An additional 
support was foreseen in front of the detector to install 
a 50 cm long steel shadow bar designed for background 
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measurements . 

3. Liquid scintillator detectors 

3.1. Detector description 

Neutron detectors were 4 cm long Aluminum cylindri- 
cal cells filled with liquid scintillator. In order to measure 
neutron yields at two different distances from the target, 
keeping the same solid angle, two different diameters were 
chosen. Four cells used as far detectors had the diameter 
of 4.6 cm (large detectors) while the other four, used as 
near detectors placed at half distance from the target, had 
the diameter of 2.3 cm (small detectors). The drawings of 
the two types of cells are shown in Figs. [3] and 21 





Figure 3: Drawing of large cell. 



Figure 4; Drawing of small cell. 



To enhance light collection, the cell Aluminum walls 
were painted with white diffusive paint EJ520 (analog of 
NE-561 and BC-622A). SmaU ceUs with Teflon diffuser 
were also produced for comparison. Reflectivity of the 
paint varies from 83% at 400 nm to 93% at 520 nm. Teflon 
was suggested to have somewhat better reflectivity (l6j . 
however we did not observe any significant difference. 

The cells were filled in Ar/N2 atmosphere with EJ301 
liquid scintillator, having a good PSD capability. The scin- 
tillator has density p = 0.874 g/cm^, atomic H:C ratio 
1.212, refractive index no = 1.505, it produces 12,000 blue 
photons per 1 MeV electron equivalent (MeVee) deposited 
energy with the first three excitation components decaying 
in 3.16, 32.3 and 270 ns. 

The liquid scintillator did not fill the entire cell volume, 
leaving a small bubble of Ar/N2 for a possible thermal 
expansion. The total amount of scintillator in each cell was 
measured by a comparison of full to empty cell weights. 



The measured expansion volume was compatible with the 
designed value of 4%. 

Each cell was sealed by a 4 mm thick borosilicate glass, 
coupled to light sensor by optical grease. As light sen- 
sors we selected Electron Tubes Enterprise PMTs 9954 
with high linearity voltage dividers C649. The choice was 
mainly determined by a good overall quantum efficiency, 
extended in the red region. This was particularly impor- 
tant for PSD 
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3.2. Calibrations 

In order to calibrate the signal of the liquid scintilla- 
tors in terms of electron equivalent deposited energy, the 
backward Compton scattering was used. This reaction al- 
lows to observe a pronounced peak of Compton scattering 
at 180 degrees and to remove underlying background by 
a fast coincidence. We used a BaF2 crystal 18), to detect 
the recoiled 7 at backward angles. The distance between 
the BaF2 crystal and the EJ301 cell was about 10-15 cm. 

The measurements were performed at LNS with a num- 
ber of 7 sources. Small detectors covering the energy range 
up to 0.5 MeVee were calibrated with ^^^Eu, ^^Na, ^"'^Bi 
and ^^''Cs sources. For the large detectors, the energy scale 
was calibrated by means of ^'^''Bi, ^•^''Cs, ^"Co and ^^Na 
sources. 

Events measured in coincidence between liquid scin- 
tillator detector and BaF2 crystal were analyzed off-line 
to select backward Compton scattering only. The TDC 
distributions of the liquid scintillator-BaF2 time-difference 
were analyzed in order to select coincidence events. Then 
we identified the corresponding peak of scattered 7 in BaF 
detector deposited energy and used the events in the peak 
for calibration. The charge distributions of liquid scintil- 
lator signals were fitted by a Gaussian to determine the 
backward Compton scattering peak position and width. 

For small detectors the 60 keVee photoelectric peak 
from an ^'*^Am source was also measured to cover the low 
energy domain. These measurements were found to be in 
agreement with coincidence data. 

Examples of obtained calibration fits for small and 
large detectors are shown in Figs. [5] and [6] The employed 
set of calibration sources allowed to extend the calibra- 
tion data almost to the threshold. The resolution of all 
detectors was found to be better than 5% at 1 MeVee. 

3.3. Efficiency evaluation 

The efficiencies of neutron detectors were obtained by 
means of Geant4.9.5-p01 Monte Carlo [l^ simulations and 
verified in a dedicated measurement. The complete detec- 
tor geometry was implemented into the simulation. Most 
of the materials were obtained from the Geant4 NIST 
database. Electromagnetic physics list was taken from 
G4EmStandardPhysics-option3 and G4EmExtraPhysics. Hadronic 
and optical physics were also activated. The scintilla- 
tor light yield, optical parameters of scintillator, reflec- 
tor painting, glass and optical grease were all configured 
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Figure 5: Small detector QDC calibration with 7 sources (upper) 
and its energy resolution (lower). 
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Figure 6: The same as in Fig.[5l but for large detectors. 



accordingly. In particular, the reflector painting was de- 
scribed as dielectric-dielectric optical surface using the uni- 
fied model and ground-front-painted finish with ctq — 10. 

A uniform neutron spectrum was generated from to 
62 MeV. The detector efficiencies were obtained as the 
ratio of detector hits depositing a visible energy above 
the experimental threshold over the number of generated 
neutrons. Since the evaluation relies on the visible en- 
ergy, a careful check of electromagnetic energy deposition 
and its saturation was performed. The comparison of the 
standard Geant4 deposited energy in liquid scintillator to 
SRIM2012 calculations [2^ showed significant disagree- 
ment for various particles. We assumed that SRIM2012 
calculations were more reliable and implemented its en- 
ergy loss tables for all relevant particles into Geant4. The 
scintillation light quenching was activated with Birks pa- 
rameter ks = 0.007 g/MeV/cm^ obtained from a com- 
bined fit of all existing world data and the data from the 
present experiment as shown in Fig. [T] The introduction 
of a second parameter leads to a much better fit quality, in 
particular in the low energy region. However, a quadratic 
term like in Ref. 21| introduces an unphysical "unquench- 
ing" behavior when extrapolated to high ionization density 



of ^^C ions. A better description of the low-energy region 
was obtained with the following parametrization: 

r^f dE \2 



AL 



1 



AE 



(1) 



with parameter C = 0.8 x lO"'' (g/MeV/cm^)^. 

The simulated efficiency has been validated in dedi- 
cated measurement, performed with the same experimen- 
tal setup used for detector calibrations. The technique 
developed in Refs. 22, was employed. This consists 
of a coincidence measurement of the neutron and gamma 
from the same fission of a ^^^Cf source. 

The obtained efficiencies for large and small detectors 
are shown in Fig. [8] and HI respectively. The measured ef- 
ficiency was found to be in good agreement with Geant4 
simulations over entire accessible energy range. This re- 
sults make us confident about using the Geant4 simulated 
efficiency in the data analysis. 

In order to avoid large systematic uncertainties only 
the energy region where the efficiency exceeded 3% was 
selected for the data analysis. 
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Figure 7: (color online) Fit of world and this experiment data on the 
light yield in the liquid scintillator. 




MeV 

Figure 8: (color online) Measured and simulated efficiency of large 
liquid scintillator detector with the threshold of 263 keVee. 




MeV 

Figure 9: (color online) The same as in Fig. |8]but for small detector 
with the threshold of 54 keVee. 



4. Data analysis 

4.I. Event selection 

At each angular setting the data were taken in few 
runs. In order to check the consistency of aU the data 
and exclude runs with significant variations in detector 
conditions, a set of quality checks were performed. These 
checks include: 

1. TDC and QDC pedestals, 

2. total number of triggers normalized to the collected 
beam charge, 

3. prompt photon peak position. 

Pedestals of TDC and QDC were stable except for some 
runs of large detectors installed at and 5 degrees. 

The total number of triggers normalized to the col- 
lected beam charge varied from 0.1% to 13% with average 
of 2.5%. We selected only those runs where the variation 
was below 5%. 

The stability of the prompt photon peak position was 
verified over each set of runs. The runs where the prompt 
photon peak position resulted to be different from the av- 
erage were discarded. 

Two independent DAQ systems were used to for small 
and large detectors. Each DAQ stored the data from four 
detectors, installed at different angles. Identification of 
trigger channel was performed during the off-line analy- 
sis. The first detector hit, corresponding to largest TDC 
value above pedestal, identified the trigger channel. The 
probability of simultaneous hits in more than one detector 
working at a few kHz overall rate within the 250 ns time 
interval covered by acquisition Gate was negligible. 

4-2. Particle identification 

Neutron hits in the detectors were identified within the 
data by a series of conditions. These included the following 
requirements: 

1. the hit time with respect to the RF signal of the 
beam had to be smaller than the RF period of 125 
ns, moreover also the amplitude of TFC signal was 
tuned to cover only the measurable range; 

2. the hit time had to be larger than the minimum time 
necessary to a neutron of 62 MeV to travel the dis- 
tance from the target to detector minus 3cr of timing 
resolution, obtained by fitting the prompt photon 
peak; 

3. the electron equivalent energy deposited in the liquid 
scintillator had to be greater than zero and smaller 
than the maximum electron equivalent energy pro- 
duced by interaction of a neutron with measured 
Time-of-Flight (kinetic energy) EXep'^iTn). 

4. deposited energy had to be greater than PSD appli- 
cability threshold. The thresholds were determined 
by the condition that both neutron loss and gamma 
contamination were below 1%. The PSD thresholds 
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were significantly higher than the hardware (discrim- 
inator) thresholds and varied from 44 to 54 keVee for 
small detectors and from 175 to 308 keVee for large 
detectors. 

5. the integrated slow scintillation component has to 
be higher than the fitted PSD separation curve. The 
determination of the optimal cut on the contribution 
of the signal tail in the total was obtained by fitting 
the PSD distributions for 7s and neutrons shown in 
Figs. [TUl and [TTl To this end the data taken without 
beam were analyzed in order to obtain 7 PSD distri- 
bution. This procedure was repeated using the data 
taken with beam on, by fitting both 7 and neutron 
distributions simultaneously. 




QDC total channel 

Figure 10: (color online) PSD distribution in large detector. The 
green dashed line shows the minimal PSD threshold; the magenta 
curves are fitted gamma peak position and width (±3cr); the blue 
curves are neutron peak position and width; the black curve repre- 
sents the cut between gammas and neutron minimizing the sum of 
inefficiency and relative contamination. 




QDC total channel 

Figure 11: (color online) The same as in Fig. 1101 but for small de- 
tectors. 



Examples of contribution of each of these conditions 
are shown in Fig. 1121 The timing cuts mainly remove the 
prompt photon peak. The deposited energy cut reduces 
the event yield at small TDC values corresponding to low 
energy neutrons. This is because background events usu- 
ally have deposited energy larger than allowed for a neu- 
tron of such small kinetic energy. PSD threshold cut leads 
to the most severe reduction of events and to the corre- 
sponding loss of efficiency, in particular in the low neutron 
energy range. PSD cut removes the remaining gamma 
background, as it can be seen from the strong reduction of 
events just before the prompt photon peak. The minimal 
efficiency cut simply limits the measurable TDC range. 
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Figure 12: (color online) Effect of different PID cuts on the TDC 
distributions measured with large detector at 30 degrees (upper) and 
small detector at 60 degrees (lower): black - all the data; red - after 
timing cuts; green - after timing and deposited energy cuts; blue - 
after timing, deposited energy and PSD threshold cuts; magenta - 
after timing, deposited energy, PSD threshold and PSD cuts; cyan 
- after timing, deposited energy, PSD threshold, PSD and minimal 
efficiency cuts. 



^.3. Collected Beam Charge Normalization 

The charge deposited by the proton beam on the Beryl- 
lium target was read out via direct electrical contact us- 
ing a high precision current digitizer. During the experi- 
ment two different current digitizers were employed: model 
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lOOOC of Brookhaven Instruments Corporation with Low 
Current Modification and Ortec model 439. Test bench 
measurements showed that the precision of these two mod- 
ules is better than 2 pA. We consider this value as a con- 
stant contribution to the normalization systematic uncer- 
tainty. Cross checks of measured charges were performed 
exploiting the linear dependence of the total event rate on 
the beam current by extrapolating the the beam current 
value to zero event rate to obtain an estimate of the cur- 
rent offset. An example of such cross check measurement 
is shown in Fig. [T21 The offset currents obtained in this 
way, ranging from -3 to pA (Ortec 439) and to 6.4 pA 
(BIC lOOOC), were used to correct the final data. The 
values of this correction were taken as the normalization 
systematic uncertainty. 




-10 10 20 30 40 50 

^beam 

Figure 13: (color online) Total event rate dependence on the mea- 
sured beam current. Different colors indicate four independent de- 
tectors installed at different angles. Straight lines show linear fit to 
the data. The extrapolated beam current at the intercept with R = 
axis (offset current) is compatible with zero. 

From the measured total beam charge deposited on 
the target Qbeam = Ibeam{i)dt the total number of 
protons on target was obtained. This was performed cor- 
recting the collected charge for the DAQ live-time fraction 
Ldaq- 

where e is the elementary charge and Ljjaq is the live- 
time correction. The latter was obtained in each run as 
the ratio between the number of recorded events and the 
number of trigger events measured by the scaler, corrected 
for multiple triggers. The measured DAQ dead time of 13 
fj,s resulted to be in agreement with expectations. 

4-4 ■ Background subtraction 

Neutrons scattered from surrounding materials, like 
beamline or detector supports, experimental hall floor or 
walls and other installations nearby, may eventually get 



into detectors. These hits will be registered as good neu- 
tron signals if the deposited energy will remain below E'^p^ . 
As mentioned before, in this experiment particular atten- 
tion was paid to reduce the material around the detec- 
tors. Nevertheless, some fraction of events comes from 
environment-scattered neutrons. In order to measure these 
environmental neutrons, special runs were taken by screen- 
ing the target direct view from the detector with a shadow 
bar. The shadow bar consisted of a 50 cm long 4.6 cm 
diameter steel cylinder. MCNP/Geant4 simulations indi- 
cated that such a shadow bar allows to reduce the initial 
neutron flux in the energy range from 0.5 to 62 MeV by 
a factor < 10"'^. Therefore, the data taken with shadow 
bar measured the contribution of environmental neutrons. 
Indeed, after subtraction of the shadow bar yields, for ex- 
ample shown in Fig. 1141 the neutron yield outside the phys- 
ical range (above prompt 7 peak) goes to zero. The con- 
tribution of this background varies from a few percent in 
average to > 50 % at the endpoints of measurable energy 
interval. 

Another background was due to PMTs dark current 
noise. This was measured during the runs without beam. 
For the large detectors this contribution was negligible. 
However, for small detectors whose discriminator thresh- 
old was set to fairly low value this background could reach 
a few percent. The contribution of this background is 
shown in Fig. [T3]in comparison with environmental back- 
ground. 

4-5. Differential neutron yield 

The differential neutron yields were extracted from the 
number of measured neutron hits in the detector A^^^*- (T„, 0„) 
according to the following equation: 

^ (T )^ ^ NffjTM 

NpdT,,dn.n^ en{Tn)Np AT„Ar!„ ' ^ ' 

where Np is the total number of incident protons on tar- 
get, en{Tn) is the neutron detection efficiency obtained in 
Geant4 simulations, T„ and 9n are neutron kinetic energy 
and production angle with respect to the beam axis, re- 
spectively. AT„ is the neutron energy bin size and Ai7„ is 
the detector solid angle. 

In order to cover the entire measured energy range we 
combined the data from the small and large detectors. The 
maximum deposited energy measurable with the small de- 
tectors was 0.5 MeVee, corresponding to maximum neu- 
tron energy of about 2 MeV. Neutrons with energies above 
2 MeV were also measured in small detectors, it was impos- 
sible to apply PSD and maximum deposited energy cuts for 
these events. On the other hand, the efhciency of the large 
detectors dropped rapidly below 2 MeV, leading to a large 
systematic uncertainty on the observed yield. Therefore, 
we selected data from the small detector at T„ < 2 MeV 
and large detector data at T„ > 2 MeV. In the overlap re- 
gion the two detector types show a fairly good agreement 
with average deviation of the order of 10%. 
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Figure 14: (color online) Shadow bar background and noise subtrac- 
tion on the TDC distributions measured with a large detector at 30 
degrees (upper) and a small detector at 90 degrees (lower): All cuts 
were applied to the data except for the timing conditions. The his- 
tograms show the following normalized spectra: black - production 
data; red - data taken with shadow bar obscuring the target view; 
green - noise measured with beam off; magenta - after background 
subtraction. 
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Figure 15: (color online) Angular dependence of the neutron yield 
for different neutron energies. 
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Figure 16: (color online) Energy dependence of the neutron yield for 
different production angles. 



The obtained yields are shown in Figs. [12] and \W\ as 
a function of neutron angle and energy, respectively. At 
low neutron energy the yield shows an isotropic behavior, 
while becoming forward-peaked at for higher energy. The 
energy dependence is well described by exponential fall-off 
with a high energy knee observed at 0„ < 60°. 

Statistical uncertainties on N^'^*- and e„ were combined 
in quadrature. The two systematic uncertainties described 
in section 14.61 instead were kept separate except for the 
fully integrated yield for which they were combined in 
quadrature. 

The comparison of the measured yields to world data 
at beam energies close to the 62 MeV of this experiment 
is shown in Fig. [TT] The data from Ref. !7|, taken at al- 
most the same beam energy, exhibit the same high energy 
structure in the measured yield at small angles. The data 
from Refs. [^, 0| taken at lower beam energies, while ex- 
hibiting similar overall behavior, are suppressed at high T„ 
due to the shrinkage of the available phase space. More- 



over, the data from Ref. (HI, which extend to lower neutron 
energies, exhibit a much milder ascent towards small T„, 
incompatible with our data. Instead the data from Ref. Q 
taken at 113 MeV have similar behavior to our points, but 
predictably lie slightly above. 



Our data were also compared to MCNP [IJ] , FLUKA [24 , 
[2^ and Geant4 [l^ simulations performed using different 
physics libraries. The results are shown in Fig.[THl MCNP 
simulations have a different descending trend as compared 
to our data such that they underestimate significantly (by 
a factor of 2 at small angles) the low energy part, while 
overestimating the high energy behavior at large angles. 
Very similar distributions were obtained with FLUKA, 
which at small angles had an overall higher magnitude. 
Geant4 provides a broad choice of physics models to test 
against our data. However, considering the energy range of 
the present experiment the set of applicable models is re- 
duced to four libraries: PreCompound model. Binary Cas- 
cade, Bertini Cascade, LEProtonlnelastic and INCLXX. 
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Figure 17: (color online) Comparison of the measure yield to the world data from Refs. 0i@>@il3- The hatched histograms show systematii 
uncertainties due to efficiency (right hatched) and normalization (left hatched). 



The PreCompound model, mostly used to describe final 
state emission, was proposed also as stand-alone model. 
Although it manages to form a sufficient number of 
compound nuclei in the p +^ Be interaction, the excess 
energy, instead of going into nucleus excitation, gets lost 
via emission of few high energy 7s. As a result the simu- 
lated yield is two order of magnitude lower than the data. 
LEProtonlnelastic does not generate two major channels: 
^Be and 2 a. It also frequently releases the excess energy 
through a high energy 7. Binary Cascade has too high 



threshold (of about 30 MeV) for ^Be production. INCLXX 
model does not have this problem, but shows similar re- 
sults. The best Geant4 model for our kinematics, Bertini 
Cascade, gives results similar to MCNP and FLUKA, and 
also underestimates significantly the yield at low neutron 
energy. 

4-6. Systematic uncertainties 

There are two major contributions to the systematic 
uncertainties in this experiment. These were the overall 
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Figure 18: (color online) Comparison of the measured yields to the 



10 



10 



MeV 

MCNP, FLUKA and Geant4 simulations at different angles. 



normalization of the data and the knowledge of the liq- 
uid scintillator efficiency. We assumed that all other un- 
certainties were either negligible or folded into these two 
contributions. 

The normalization was based on the charge accumu- 
lated on the Beryllium target and read out by a digital 
integrator. The precision of the integrator was better than 
2 pA, but in some runs corrections to the measured charge 
were necessary as discussed in section 14.31 Therefore we 
evaluated the normalization systematic uncertainty as the 
maximum between 2 pA value and the correction applied. 

The efficiency systematic uncertainties were evaluated 
from the comparison of Geant4 simulations to the ded- 
icated measurement described in section 13.31 From this 
comparison the standard deviation of the data from Geant4 
was extracted in the range 0.5-3 MeV for the small de- 
tectors and 2-10 MeV for the large ones. The deviation 
resulted to be 10% for both types of detectors. This value 
was taken as an estimate of efficiency systematic uncer- 
tainty. 

4-7. Low energy detectors 

At low neutron energy T„ < 0.47 MeV liquid scintil- 
lators were not efficient enough to determine the neutron 



flux. This energy range was measured by two different de- 
tectors: one was a "^He tube (Canberra 150NH50/5A) and 
the other a silicon detector (Micron MSX09) covered with 
a ^Li converter. Both preamplified detector pulses were 
amplified by an Ortec 672 spectroscopy amplifiers and then 
digitized by the same DAQ system using a SILENA 4418 
peak sensing ADC. In few special runs these detectors were 
installed on the same supports used for liquid scintillators. 
Assuming the low energy neutron fiux to be isotropic we 
measured only few angles: 70 and 90 degrees for '^He tube 
and 130 degrees for ^Li-silicon detector. The event rates 
were fairly low, hence the distance to the target was cho- 
sen to be 64 cm for "^He tube and 69 cm for ^LiF-silicon 
detector. 

The ■^He tube energy response to thermalized neutrons 
from an AmBe source was first analyzed for the purpose 
of ADC calibration. Using this calibration, energy spectra 
of the data taken during the experiment were obtained. 
In order to remove environmental background of neutrons 
scattered from the surrounding materials shadow bar mea- 
surements were subtracted from the data. Since the "^He 
tube geometry was not optimized for this particular ex- 
periment, the contribution of environmental background 
resulted to be very significant leading to large systematic 
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uncertainties. The obtained spectrum was divided in five 
bins from to 0.47 MeV according to the energy resolu- 
tion estimated from Geant4 simulations. Then an unfold- 
ing procedure was applied to the data. First we simulated 
in Geant4 the ^He tube response to the neutron flux mea- 
sured with liquid scintillators and subtracted it from the 
data. Then we extrapolated this flux to lower energies 
via dependence, with f3 being an adjustable parame- 
ter, and extracted neutron yield in defined kinematic bins 
using Geant4 simulations of efficiency. This step was re- 
peated few times adapting at every iteration the simulated 
neutron yield to the extracted one. Once the extraction 
procedure converged, we obtained the final neutron yield 
at low energies shown in Fig. [TO] in comparison to liquid 
scintillator data. 
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Figure 19: (color online) Energy dependence of the neutron yield 
at low energy measured by the "^He tube in comparison with liquid 
scintillator data. Red triangles represent •^He data, while red dashed 
area gives its systematic uncertainty. Black squares and green circles 
show liquid scintillator data with and without PSD. 



The highest point of ^He data was found to be in good 
agreement with liquid scintillator measurements. In par- 
ticular, the data from liquid scintillator without PSD ex- 
tend down to 0.3 MeV with somewhat larger systematic 
uncertainties. Although the estimated systematic uncer- 
tainties of ^He data resulted to be > 50%, the obtained 
data allowed to constrain the missing part of the neutron 
energy spectra. The systematic uncertainty of ^He data 
contributes for 50% to the uncertainty of the integrated 
neutron yield, although in absolute value it corresponds 
to only 3% of the yield. 

The ^LiF-silicon detector data confirmed the neutron 
yield extracted from liquid scintillator measurements and 
put compatible upper limits on '^He data. However, due to 
large silicon active layer thickness it was mainly sensitive 
to high energy neutrons. 

4-8. Integrated yield 

Combining liquid scintillator data and '^He tube data 
the complete neutron energy coverage was achieved. Thus 



the integrated neutron yield was calculated separately for 
each angular setting. This was performed assuming that 
the low-energy part of the spectrum was isotropic, sup- 
ported by the low energy points from liquid scintillator 
detectors. The obtained angular distribution of the inte- 
grated neutron yield is shown in Fig. 1201 
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Figure 20: Angular dependence of the integrated neutron yield. The 
hatched histograms show systematic uncertainties due to efficiency 
(right hatched) and normalization (left hatched). 



The integrated neutron yield falls rapidly with the an- 
gle up to about 90 degrees, while at larger angles the dis- 
tribution is almost flat. This suggests a large contribu- 
tion of the process where the incident proton kicks off the 
loosely bound neutron via a quasi-elastic scattering. In- 
stead, at large angles the dominant process is, perhaps, 
of evaporative nature. By integrating, in turn, the neu- 
tron yield in Fig. EUl over the angle, the total neutron yield 
Nn/Np = 0.1003 ± 0.0003,tat. ± 0.0053,.ys. n/p was ob- 
tained. This value can be compared to the MCNP simula- 
tions iV„/7Vp = 0.103 n/p (ENDF VII) or N^/Np = 0.096 
n/p (LA150) performed for the same beam energy, and to 
the experimental data from Ref. [i| 7V„/7Vp = 0.110±0.007 
n/p at 65 MeV. The total yield from Ref. ijil was measured 
by thermalization of produced neutrons and their subse- 
quent capture. 

5. Conclusions 

Fully differential neutron yields produced by 62 MeV 
proton beam impinging on a thick, fully absorbing ^Be tar- 
get were measured in almost complete energy and angular 
ranges. The precision of the data is limited by systematic 
uncertainties on the detector efficiency and by the abso- 
lute normalization. These uncertainties of the yields are of 
the order of 10%. A better precision on detector efficiency 
could be achieved by means of a precisely calibrated neu- 
tron source, while the normalization uncertainty could be 
improved by using a higher beam intensity. 

The obtained data are in agreement with previous ex- 
periments in 60-70 MeV beam energy range, extending 
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considerably the kinematic coverage. The overall yield, in- 
tegrated over neutron energy and angle, is in good agree- 
ment with a measurement based on thermalization and 
capture of produced neutrons Q • The comparison of our 
data to MCNP and FLUKA simulations showed some de- 
viations in angular and neutron energy behaviors, although 
the the integrated yield values were found to be consistent. 
The comparison to various Geant4 models exhibited their 
limitations in the low beam energy domain. 

The presented data will be used in the design of the 
ADS core Q to describe its neutron source. These data 
can also be useful for the development of the neutrino 
source proposed recently in Ref. [26|. 
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